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ABSTRACT 

We present the results of Hubble Space Telescope Wide Field Camera 3 observations of the core of 
the Phoenix Cluster (SPT-CLJ2344-4243) in five broadband filters spanning rest-frame 1000-5500A. 
These observations reveal complex, filamentary blue emission, extending for >40 kpc from the bright- 
est cluster galaxy. We observe an underlying, diffuse population of old stars, following an r 1 / 4 distri- 
bution, confirming that this system is somewhat relaxed. The spectral energy distribution in the inner 
part of the galaxy, as well as along the extended filaments, is a smooth continuum and is consistent 
with that of a star-forming galaxy, suggesting that the extended, filamentary emission is not due to 
a large-scale highly- ionized outflow from the central AGN, but rather a massive population of young 
stars. We estimate an extinction-corrected star formation rate of 798 ± 42 M yr -1 , consistent with 
our earlier work based on low spatial resolution ultraviolet, optical, and infrared imaging. We argue 
that such a high star formation rate is not the result of a merger, as it would require >10 mergers with 
gas-rich galaxies and there is no evidence for such multiple merger events. Instead, we propose that 
the high X-ray cooling rate of ~2850 Mq yr -1 is the origin of the cold gas reservoir. The combination 
of such a high cooling rate and the relatively weak radio source in the cluster core suggests that 
feedback has been unable to halt runaway cooling in this system, leading to this tremendous burst of 
star formation. 



1. INTRODUCTION 

The cores of galaxy clusters represent one of the few 
places in the Universe where large-scale cooling and feed- 
back processes can be readily observed. Unlike isolated 
massive galaxies, where some of the energy injected into 
the interstellar medium (ISM) from the active galactic 
nucleus (AGN) often escapes into the low-density inter- 
galactic medium (IGM), the denser intracluster medium 
(ICM) in cluster cores retains an i mprint of this feed- 
back in the form of cav ities (e.g., Birz an et al.l 2004; 
Hlavacek-Larrondo et al. 2012) or ripples (e.g., Perseus 
A; iFabian et al.ll2003l ). Similarly, the accretion of hot 
gas from the IGM onto massive galaxies is challenging 
to observe due to the very low densities, while such phe- 
nomena have been studied in depth for decades in galaxy 
clusters. These so-called "cooling flows" in galaxy clus- 
ters, which were once thought to be massive flows of cool 
gas on the o rder of 100-1000 M yr -1 (see review by 



Fabian 1994), are now understood to be considerably less 
massive, depositing on the order of 1-10 M yr -1 of cool 
gas onto the brightest cluster galaxy (BCG). This resid- 

1 Kavli Institute for Astrophysics and Space Research, MIT, 
Cambridge, MA 02139, USA 

2 Kavli Institute for Cosmological Physics, University of 
Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 

3 Department of Astronomy, University of Maryland, College 
Park, MD 20742, USA 

4 Joint Space-Science Institute, University of Maryland, Col- 
lege Park, MD 20742, USA 

5 Astroparticle Physics Laboratory, NASA Goddard Space 
Flight Center, Greenbelt, MD 20771, USA 

6 Max- Planck- Inst it ut fur extraterrestrische Physik, Postfach 
1312, D-85741 Garching, Germany 

7 Department of Physics, University of California, Berkeley, 
CA 94720, USA 

Email: mcdonald@space.mit.edu 
t Hubble Fellow 



ual c ooling has been dete c ted in UV spectra (O VI lines; 
e.g., lOegerle et aTHgOOlUBregman et al. 2006h. warm 



ionized gas (e.g.. iHu et al.l IT9 85: Crawford et all 119991 : 
[McD onald et al. 2010, 20 lTaF wa rm m olecular gas (e.g., 
Uaffe et a l. 2005; Donahue et al. 2011 ), cold molecular 
gas (e .g^Edge 2001; Edge et al. 2002; S alome fc Combes! 
2003; McDo nalcTet al.l 120 12bF7^ and via young star s de- 
tected in the far UV (e.g. , McDonald et al. 2011b) and 
mid infrared (e.g.. lO'Dea et al . 2008; Rawle et al. 2012h. 
The exact feedback mechanism that prevents the cooling 
catastrophe from happening, allowing only a few per- 
cent of the predicted cooling flow to reach low (< 10 4 K) 
temperatures, is still not fully understood. The leading 
hypothesis of mechanical AGN-driven feedback is sup- 
ported by the strong correlation between the total energy 
required to offset cooling and the amount of radio power 
required to inflate X-ray cavities (McNamara fc Nulsenl 
120071 : lFabiaiJl2Q12t iMcNamara fc N ulsen 2Q12f ). 

Recently. iMcDonald et al.l (j2012af ) reported the unique 
properties of a galaxy cluster at z = 0.596, SPT- 
CLJ2344-4243 (hereafter the Phoenix Cluster), which 
was initially discovered by the Sout h Pole Telescope 
using the Sunyaev Zel'dovich effect (Wi lliamson et al.l 
12011 . This cluster is among the most massive 
(M 2 oo,y x ~ 2 -5 x 10 15 M ) and X-ray luminous 
(L2-iokeV = 8-2 x 10 45 erg s _1 ) yet discovered. In the 
central ~100 kpc, there is a dense, highly-luminous cool 
core, with a classical cooling rate (^ = 2 L £j^ v ) higher 
than any known cluster, at 3820 M yr -1 . Surpris- 
ingly, while most nearby cool core clusters convert only 
a few percent of the classical cooling flow into stars, the 
central galaxy in the Phoenix cluster appears to be ex- 
periencing a 740 M yr -1 starburst, corresponding to 
roughly 20% of the predicted cooling flow. However, 
the central galaxy also contains a powerful AGN which 
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Fig. 1. — HST-WFC3 images of the core of the Phoenix cluster, in 5 different optical bands (rest-frame 1000-5500A). These images 
highlight the complex morphology of the BCG. The presence of extended emission in all bands argues in favor of continuum, rather than 
line, emission. The full extent of these filaments, ~40 kpc, is reminiscent of those observed in nearby cool core clusters such as Perseus 
and Abell 1795. 



makes it difficult to separate contributions to the ultra- 
violet, optical, and infrared emission from star- forming 
regions versus the AGN in low-resolution imaging. Fur- 
ther, the available ground-based data optical and space- 
ba sed ultraviolet (UV) and infrared (IR) data presented 
in iMcDonald et al.l (j2012a[ ) were unable to distinguish 
between in situ star formation (i.e., cooling flows), or 
the late stages of a merger. 

In this Letter, we present new Hubble Space Tele- 
scope observations which improve significantly in depth 
and spatial re s olution on the data presented in 
IMcDonald et al.l (|2012afl . In §2 we describe these data 
and our analysis, after which we present our new, de- 
tailed view of the Phoenix cluster in §3. In §4 we dis- 
cuss the possible interpretations of these data, including 
whether or not we can separate the cooling flow scenario 
from a pure AGN or merger scenario. We conclude with 
a summary of these results in §5. Throughout this letter 
we assume Ho = 70.2 km s _1 Mpc -1 , VLm = 0.278, and 
Q A = 0.728. 

2. HST DATA 

To study in detail the purported starburst in the core 
of the Phoenix cluster, we obtained broadband imaging 
with the Hubble Space Telescope Wide Field Camera 3 
(HST WFC3) in five optical filters - F225W, F336W, 
F475W, F625W, F814W - which span rest-frame wave- 
lengths from ~1000A to ~5500A. These observations 
were carried out over 2 orbits of Director's Discretionary 



Time, with an average exposure time of ~ 1800s per filter 
(PID #13102, PI McDonald). 

In each filter, a series of two dithered exposures were 
obtained. The LA CosmicE3 (|van Dokkuml [2QQlh soft- 
ware was run on individual exposures, generating ac- 
curate cosmic ray and bad pixel masks for each image. 
Pairs of exposures were then combined using the PyRAF 
astrodrizzle routin^E], with the aforementioned masks 
preceding the standard cosmic ray detection in the MUL- 
tidrizzle task. The final, cleaned images are presented 
in Figure [TJ 

All optical and U V fluxes were c orrect ed for Galactic 
extinction following iCardelli et al.l (|1989f ) using a Galac- 
tic reddening estim ate of E(B — V ) = 0.0 17 towards the 
cluster center, from Schleg el et al.l (|1998[ ). 

3. RESULTS 

In Figure [1] we show the newly- acquired far-UV 
through optical HST images of the core of the Phoenix 
cluster, which provide a detailed picture of this ex- 
otic system. These images show significant, extended 
filamentary emission at all wavelengths from ~1000- 
5500A, overlaid on a relatively smooth population of 
older, red stars. The most extended pair of filaments to 
the north of the BCG are (40 kpc) in length, sim- 
ilar to the most extended filaments seen in Abell 1795 



http : //www . astro . yale . edu/dokkum/lacosmic/ 

http : //www . stsci . edu/hst/HST_overview/ drizzlepac 



Fig. 2. — Color image, combining the F475W, F625, and F814W bands, showing the young, filamentary, star-forming regions overlaid on 
the diffuse, old stellar component of the BCG. In the inset we show the central ~20 kpc, with the position (and uncertainty) of the X-ray 
point source in green. There appear to be gaps in the emission to the northwest and southeast of the peak, possibly due to strong dust 
lanes. 



( McD onald & Veilleux 2009), and the Pe rseus cluster 
( Conselice et all 120011 : iFabian et all [2008). We mea- 
sure a total rest-frame far-UV flux density of fF225W 
= 1.26xl0 -27 erg s _1 cm -2 Hz -1 , consistent with 
the GA LEX-derived flux presented in iMcDonald et al.l 
([2012^ . 

The fact that such complex, filamentary morphology is 
present in all five filters suggests that the BCG is forming 
stars at a prodigious rate. In the wavelength range cov- 
ered, there may be contributing emission from the Lya 
(F225W), [O II] (F625W), and [O III] (F814W) lines. 
However, the F336W and F475W bands, which have sim- 



ilar surface brightnesses to the other three bands, should 
be relatively free from emission lines, suggesting that ion- 
ized gas is not the dominant source of the observed flux 
in Figure [TJ 

In Figure [2] we show a three-color (F475W, F625W, 
F814W) image of the cluster core. This figure shows a 
clear difference in the stellar populations between the 
young (blue) filaments and the underlying, smoothly- 
distributed, old (red) stars. The peak of the emission in 
all bands is coincident (within the positional uncertain- 
ties) with the X-ray point source. To the northwest and 
southeast of the emission peak are dark lanes, most likely 
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Fig. 3. Left: Surface brightness profiles of the BCG in all five bands (purple: F225W, blue: F336W, green: F475W, orange: F625W, 
red: F814W) . Each profile represents the average along four radial cuts with orientations chosen to avoid most of the extended, filamentary 
emission. The diffuse, red component of the BCG has a surface brightness distribution typical of an elliptical galaxy, well model ed by a r 1 / 4 
fall-off in surface brightness (dashed line; de Vaucouleurs 1948). Middle: 2-dimensional surface brightness model from GALFIT ( Peng et al. 
2010), with ellipticity and PA of 0.88 and -52°, respectively. Right: Residual image generated by subtracting the r 1 / 4 model from each of 
the three reddest bands. Both color images use the same color scaling as Figure [J] This figure highlights the complex morphology of the 
star-forming filaments superimposed on the smooth, old stellar population in the BCG. 



due to obscuration by dust. Overall, the color of the fila- 
mentary emission appears roughly constant with radius, 
and is reminiscent of a young, star-forming galaxy. We 
see no evidence for multiple bulges or tidal features, both 
of which would indicate that this system is the result of 
a recent merger of gas-rich galaxies. 

In addition to the complex, filamentary emission, we 
observe two radial features to the northwest of the cen- 
tral galaxy, at a location >60 kpc from the X-ray peak. 
We interpret these as two radial arc candidates, back- 
ground galaxies gravitationally lensed by the BCG of the 
Phoenix cluster. A more thorough analysis of the cluster 
potential is required to assess this interpretation, which 
is beyond the scope of this letter. 

Figure [3] we shows the multi-band surface brightness 
profiles of the BCG (left panel), which have been com- 
puted along radial cuts at four different angles (90°, 120°, 
180°, 210°), chosen to avoid the blue, filamentary emis- 
sion. The profile follows a r 1 / 4 profile, which is typical of 
relaxed, early- type galaxies ([de Vaucouleurs|[l948f ). Pro- 
jecting these 1-dimensional profiles back onto the sky, 
we can separate diffuse, giant elliptical emission (middle 
panel) from clumpy, star- forming emission (right panel). 
This confirms that the star-forming filaments do not re- 
semble merging disk galaxies: there is a notable absence 
of smooth, arcing tidal features and multiple-bulges. We 
note that such r 1 / 4 surface brightness distributions are 
also common in the fi nal stages (single-nucleus) of gas- 
rich mergers fULIRGs: fVeilleux et alJl2002L I2006D . How- 
ever, in contrast to local ULIRGs, this BCG is a fac- 
tor of ~4 larger in size (IVeilleux et al.ll2002f), a factor 
of ^60 higher in stellar mass (|Veilleux et all 12006( 1. and 
resid es in an environment ^50-100 times richer than nor- 
mal ( Zauderer et al.l l2007). These factors argue that the 
Phoenix BCG is unlike a traditional local ULIRG by any 
definition other than the high total infrared luminosity. 

In Figure H] we provide the spectral energy distribution 
(SED) in several representative regions around the BCG. 
The diffuse emission at large and small radii indicate a 
significant positive age gradient in the diffuse population. 



At large (r ~ 40 kpc) radii, the diffuse emission is consis- 
tent with a 2-5 Gyr old elliptical galaxy, while at smaller 
radii (~20 kpc) the diffuse stellar populations are much 
younger, similar in color to a star- forming late- type spi- 
ral galaxy. The extended, morphologically-complex fila- 
ments, after subtraction of the diffuse stellar component, 
show an excess of UV emission at all radii. The SED of 
the brightest filaments appear remarkably similar to the 
diffuse component in the central region, suggesting that 
these stars are being mixed on short timescales. In the 
faintest, most extended filaments, there is a substantial 
excess of emission in the F225W filter, at the location of 
the redshifted Lya line, suggesting that these filaments 
may also contain warm (> 10 4 K), ionized gas. The over- 
all flatness of the SED in the UV-bright regions is ex- 
actly what one would expect for a mix of young stars 
and warm, ionized gas, given the width of the broad- 
band filters. At the peak of the optical emission, coinci- 
dent with the X-ray source, the SEP is w ell matched by 
a dusty type-2 QSO (Poll etta et al.l 120071 ). which is con- 
sistent with our X-ray observations of a highly-reddened 
hard X-ray point source. Finally, we find that the lin- 
ear features to the northwest of the BCG (Figure [2]) is 
inconsistent with synchrotron emission, suggesting that 
it is either an exceptionally linear filament or a gravita- 
tional arc. 

The combination of Figures [THU paint a picture of an 
old, giant elliptical galaxy that is experiencing a resur- 
gence of star formation. Below we re-evaluate the SFR 
in this system and describe various scenarios to explain 
this star-forming activit y, building on the discussion of 
iMcDonald et al.l (|2012at ). 

4. DISCUSSION 

The deep, high spatial resolution HST UV and optical 
imaging presented in §3 have revealed an exceptionally 
complex system. Below we utilize this improved spatial 
resolution to estimate a new, UV-derived SFR for the 
BCG in the Phoenix cluster, followed by a discussion of 
three possible origins for the extended, filamentary UV 
emission. 
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Fig. 4. — Left: F814W image of the Phoenix BCG, with positions of nucleus, filaments, and diffuse emission indicated. Right: SEDs for 
the different regions indicated on the left. For comparison, we also show the SEDs for four nearby early-type satellite galaxies. Orange 
and red spectra correspond to passive stellar populations of ages 2 and 5 Gyr, respectively, wh ile blue and purple spectra correspond to a 
typical late-type spiral galaxy (Sdm) and dusty type-2 QSO, respectively (Polletta et al. 2007). All model spectra are normalized by eye. 
The faintest (labels 1 and 2) filaments show evidence for nebular emission, with peaks in the F225W and F625W bands, consistent with 
Lya and [O II] emission. This figure confirms that, at large radii the stellar populations are typical for an early-type galaxy, while at small 
radii both the diffuse and filamentary emission result from significantly younger stellar populations. We also see that the UV-optical SED 
of the nucleus is consistent with a dusty AGN, while the linear feature to the northwest is inconsistent with synchrotron emission from a 
jet. 
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(|2012aD from Balmer line ratios (E(B-V) = 0.34). Using 
the local reddening correction and uncertainty we can ac- 
curately correct the UV luminosity, resulting in a more 
confident estimate of L2000 = 5.7 ± 0.3 x 1 3Q erg s -1 
and a UV-derived SFR of 798 ± 42 M yr" 1 (jKennicuttl 
fl998h . 

We note that this estimat e is consistent with the AG N- 
subtracted SFR quoted in iMcDonald et al.l (|2012af ) of 
739 ± 160 M yr -1 and with the empirical, ex tinction- 
implicit, method of iRosa- Gonzalez et al.l (|2QQ2f ). which 
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Fig. 5. Left: Image in the F475W band, with 0.2" x 0.2" bins 
overlaid in regions with > 5cr detections in both the F336W and 
F475W bands. Right: Smoothed reddening map, derived assum - 
ing a flat UV SED in the absence of reddening (Kennicutt 1998). 
This map allows us to perform local reddening corrections on the 
UV images, yielding a more tightly constrained estimate of the 
extinction-corrected UV luminosity. 

4.1. A Revised Estimate of the SFR 

The SFR reported in IMcDonald et al.l (|2012afl . while 
utilizing an array of multi-wavelength data, necessarily 
required multiple assumptions to remove AGN contam- 
ination. With the addition of high-resolution HST UV 
imaging, we avoid such assumptions, leaving only the 
(typical) assumptions of dust extinctions and star for- 
mation laws. 

To estimate the reddening, we assume that the 
reddening-free SED is flat at 1500-3000 A (jKennicuttl 
I1998D . Since the F225W filter is likely contaminated 
by Lya emission, we use the mean fF336w/fF475W (rest 
frame /210//298) flux ratio to derive a reddening correc- 
tion, assuming fF33ew/fF475W = 1 for E(B-V) = 0. In 
Figure [5] we show the spatial distribution of the redden- 
ing, E(B-V), over regions with significant (S/N > 5) UV 
flux. We note that the mean reddening i n this map agrees 
well with the reddening presented in IMcDonald et al.l 



4.2. Possible Sources of Extended UV Emission 

4.2.1. AGN-driven Outflow 

Inspired by the high far-infrared luminosity, coupled 
with the hard X-ray and radio point source, a viable 
explana tion for the data presented by IMcDonald et al.l 
(|2012af ) is a dusty AG N driving a large-scale outflow 
(e.g., IRAS 09104+4109: lO'Sullivan eTatl 120121 ) . While 
there certainly is a powerful AGN in the core of the 
Phoenix cluster, the new HST data suggest that the 
substantial UV luminosity can not be fully attributed 
to this AGN. Figure [1] shows that only a small fraction 
(<10%) of the total UV luminosity originates from a cen- 
tral point source, with the majority originating from ex- 
tended, filamentary regions. Furthermore, the high rela- 
tive fluxes in the F336W and F475W bands, which should 
be free from line emission, suggest that the majority of 
the UV/optical flux in these complex filaments is contin- 
uum emission. Finally, while the central galaxy does not 
appear to host a disk, it does appear elongated in the 
diffuse component (b/a ~ 0.8) towards the northwest 
(PA ~ -40°), in the same direction as the extended fila- 
ments. There is a general lack of UV/blue emission along 
the minor axis of the central galaxy, contrary to what 
is typically observed in wide-angle AGN-driven outflows 
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(jVeilleux et aLll2QQ5h . Thus, we argue that the filamen- 
tary UV emission does not result from an outflow. 

4.3. Gas Rich Merger (s) 

The far-IR luminos ity (L FIR = 9.5 ± 1.1 xlO 12 L ; 
iMcDonald et al.ll2012af ) of the BCG in the Phoenix clus- 
ter is significantly higher than the threshold (>10 12 L ) 
to be called an ultr alum inous infrared galaxy (ULIRG; 
iSanders fc Mirabel 1 19961 ), and pushing into "hyper lumi- 
nous" territory (HyLIRG; >10 13 L ). The majority 
of ULIRGs appe ar to be late-stage mergers of gas-rich 
galaxies f~95%: IVeilleux et al.ll2QQ2l ). 

In a merger scenario, all of the cold gas would have 
to come from one or more infalling gas-rich galaxies, as- 
suming that the giant elliptical galaxy was gas-poor to 
begin with. Following iHopkins et al.l (|2008[ ). we assume 
that the peak star formation happens during the "coa- 
lescence/ULIRG" stage, and that this star formation is 
roughly constant over 100 Myr. Assuming SFR = 800 
M yr _1 , and 100% efficiency of star formation (being 
conservative), we estimate that the total mass of cold 
gas re quired would be 8xl0 10 Mq. From iBoselli et aTl 
(1997), gas-rich spiral galaxies in the cluster environ- 
ment have typical cold gas masses of M# 2 = io 8 9±0 - 4 
Mq, implying that we would need ~100 "typical" disk 
galaxies, or ~10 exceptionally gas-rich galaxies to fulfill 
our conservative requirements for fuel. In the latter case, 
the question of how such gas rich galaxies make into the 
dense cluster core without being stripped further com- 
pounds the problem. 

These arguments, combined with the fact that this 
BCG is unlike "typical" ULIRGs in size, stellar mass, 
and environment (§3) suggest that the high SFR is not 
a result of gas-rich mergers. 

4.4. Runaway Cooling of the Intracluster Medium 

Our preferred explanation in IMcDonald et al.l (|2012a| ) 
was that the star formation in the cluster core is being 
fueled by gas cooling out of the ICM. This remains the 
most plausible avenue for such a large amount of cold gas 
to reach the core of the cluster, and is supported by the 
exceptionally bright X-ray cool core and relatively weak 
radio source. Such an imbalance between cooling and 
feedback could lead to runaway cooling, fueling bursts of 
star formation. 

Following Wh ite et al.l (jl997[ h we estimate from the 
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X-ray data the mass deposition rate of the cooling flow, 
combining the X-ray cooling luminosity with the gravita- 
tional potential of the cluster in order to correct for grav- 
itational work done as the gas falls towards the cluster 
center. We obtain an ICM cooling rate of M = 2850 M 
yr _1 which is enough to fuel a 798 M yr _1 starburst, 
assuming the feedback mechanism that prevents star for- 
mation in nearby clusters is operating less efficiently 
in Phoenix. In the future, hi gh resolution X-ray spec- 
troscopy of the cool core (e.g. JPeterson fc Fabianl f2006) 
will provide firm estimates of the ICM cooling rate down 
to low temperatures, revealing whether or not this sce- 
nario is a plausible fuel source for this massive starburst. 

5. SUMMARY AND CONCLUSIONS 

We report new HST observations of the Phoenix clus- 
ter (SPT-CLJ2344-4243) with WFC3-UVIS in five fil- 
ters covering rest-frame wavelengths 1000-5500A. The 
high spatial resolution of HST is able to separate bright 
UV emission from the AGN and the surrounding dif- 
fuse, extended emission, definitively confirming the pres- 
ence of a starburst in the BCG. The morphology of 
this central galaxy is complex, with narrow filaments ex- 
tending for >40 kpc, reminiscent of the nearby Perseus 
and Abell 1795 clusters. We argue that the majority 
of the observed UV emission is due to young stars, on 
the basis of the complex morphology and flat SED over 
the wavelength range 1000-5500A. We confirm the high 
SFR presented in Mc Donald et al.l (j2012af h measuring an 
extinction-corrected, AGN-removed UV-derived SFR of 
798 ± 42 M© yr _1 . We find that merger driven scenarios 
would require an unreasonably large number of gas rich 
galaxies to supply the cold gas reservoir required to fuel 
the starburst, and conclude that the starburst is likely 
fueled by a runaway cooling flow. 
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